Stability of Barycentric Interpolation Formulas for Extrapolation
نویسندگان
چکیده
The barycentric interpolation formula defines a stable algorithm for evaluation at points in [−1, 1] of polynomial interpolants through data on Chebyshev grids. Here it is shown that for evaluation at points in the complex plane outside [−1, 1], the algorithm becomes unstable and should be replaced by the alternative modified Lagrange or “first barycentric” formula dating to Jacobi in 1825. This difference in stability confirms the theory published by N. J. Higham in 2004 [IMA J. Numer. Anal., 24 (2004), pp. 547–556] and has practical consequences for computation with rational functions.
منابع مشابه
Stability of Barycentric Interpolation Formulas
The barycentric interpolation formula defines a stable algorithm for evaluation at points in [−1, 1] of polynomial interpolants through data on Chebyshev grids. Here it is shown that for evaluation at points in the complex plane outside [−1, 1], the algorithm becomes unstable and should be replaced by the alternative modified Lagrange or “first barycentric” formula dating to Jacobi in 1825. Thi...
متن کاملOn the backward stability of the second barycentric formula for interpolation
We present a new stability analysis for the second barycentric formula, showing that this formula is backward stable when the relevant Lebesgue constant is small.
متن کاملThe numerical stability of barycentric Lagrange interpolation
The Lagrange representation of the interpolating polynomial can be rewritten in two more computationally attractive forms: a modified Lagrange form and a barycentric form. We give an error analysis of the evaluation of the interpolating polynomial using these two forms. The modified Lagrange formula is shown to be backward stable. The barycentric formula has a less favourable error analysis, bu...
متن کاملOn the numerical stability of the second barycentric formula for trigonometric interpolation in shifted equispaced points
We consider the numerical stability of the second barycentric formula for evaluation at points in [0, 2π ] of trigonometric interpolants in an odd number of equispaced points in that interval. We show that, contrary to the prevailing view, which claims that this formula is always stable, it actually possesses a subtle instability that seems not to have been noticed before. This instability can ...
متن کاملFirst applications of a formula for the error of finite sinc interpolation
In former articles we have given a formula for the error committed when interpolating a several times differentiable function by the sinc interpolant on a fixed finite interval. In the present workwe demonstrate the relevance of the formula through several applications: correction of the interpolant through the insertion of derivatives to increase its order of convergence, improvement of the ba...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- SIAM J. Scientific Computing
دوره 34 شماره
صفحات -
تاریخ انتشار 2012